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Stanford, California 


ABSTRACT 

The purpose of this study is to compare predictions of the back- 
scatter spectrum, including effects of ionospheric inhomogeneity, with 
experimental observations of incoherent backscatter from an artificially 
heated region. Our calculations show that the strongest backscatter 
echo received is not from the reflection level, but from a region some 
distance below (about 9OO m for an experiment carried out at Arecibo) . 

By taking the standing wave pattern of the pump into account properly, 
the present theory explains certain asymmetrical features of the up- 
shifted and down-shifted plasma lines in the backscatter spectrum, and 
the several satellite peaks typically accompanying them. 
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1, Introduction 


From experimental results obtained at Arecibo ( Carlson^ Gordon and 
Showen, 1972; Kantor, 197^)^ it is believed that the enhanced heating 
and anomalous absorption occurring in the ionosphere are mainly due to 
parametric decay instabilities: the HF ordinary (pump) wave supplied by 

a high power transmitter 100 kw)^ and propagating in the F-layer^ 
decays into electron plasma and ion acoustic waves to enhance the plasma 
fluctuations ( puBois and Goldman, I965 ; Nishjkawa, I968) . Several 
analyses have been carried out of the saturation spectra in homogeneous 
( PuBois and Goldman, 1972 a^b ; ValeOj Oberman and Perkins, 1972 ; Kruer 
and Valeo, 1973; Harker, 1972; Kuo and Fejer, I 972 ; Fejer and Kuo, 1973; 
Perkins, Oberman and Valeo, 1974) and inhomogeneous ( Arnush, Fried and 
Kennel, 197^) plasmas^ with the assumption that the pump field is well 
above the threshold for parametric amplification. This results in 
saturation of the electron plasma waves due to nonlinear Landau damping 
by positive ions, or due to absorption caused by electrons whose orbits 
are nonlinearly perturbed by the waves ( Bezzerides and Weinstock, 1972). 

Some essential features of the backscatter spectrum are not fully 
understood yet, however. For example, Carlson et al . (I972) have 
obtained a pair of plasma lines in the spectrum enhanced by the pump 
wave (Iq - 5 •62 MHz). These plasma lines are shifted up and down from 
the incoherent scatter diagnostic beam ^ 430 MHz) by approximately 
~ ) where f^(?«4 kHz) is the frequency of the ion acoustic waves 

generated, and differ from each other in the intensities of the peaks 
(Figure l) . Kantor (1974) has further observed that this frequency 
displacement varies in time, but its difference for the up-shifted and 
the down-shifted lines is almost constant throughout the day, averaging 
0.4 ± 0.1 kHz (Figure 2). As shown in Figure 1, he also observed 
several maxima and minima near the peaks of the spectrum. The purpose 
of this paper is to present a possible explanation of these features by 
considering the scattering region in relation to the Airy function form 
of the pump field variation in an ionosphere with an approximately 
linear density profile. 
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TEMPERATURE (10 “K) 



FIG, 1. Backscatter spectrum obtained at Arecibo, (Af is frequency 
displacement of up-shifted line from f, + f^ , and of down-shifted 
line from f^ - f^; ~ ^30 MHz; ordinary mode pump wave^ 

f^ = 5,62 MHz. After Kantor, 197i^-) , 
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FIG. 2. Diurnal variations of frequency displacements. 

(The figure contains data obtained between May 1971 
March 1972. The bars represent the range of all available 
data for one hour intervals. After Kantor, 197^)* 
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In regard to the asymmetry of the intensities, Fejer (1972) suggested 
that it is due to the convective amplification of upward- and downward- . 
propagating plasma waves: those responsible for scattering the up-shifted 

and the down-shifted lines have been generated below and above the 
scattering height, respectively. This generation height difference was 
found to be nearly equal to the adjacent maximum -minimum distance of the 
pump amplitude. Following the estimate of Perkins et al . (19T^), however, 
we neglect this convective amplification. Instead, we consider the small 
difference between the wave-numbers of the up-shifted and down-shifted 
lines, from which it is found that the asymmetry of the line intensities 
is due to the height difference of the scattering regions for the two 
lines: the up-shifted line is scattered in a region below that for the 

down-shifted, by about 1/10 of the distance of the adjacent maximum- 
minimum of the pump amplitude. It is further shown that the present 
analysis is consistent with Kantor’s observations. 
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2 . Incoherent Scattering Regions 

As the pump wave from the HF transmitter propagates vertically 
upward in the ionosphere, its group velocity decreases very sharply 
near, and reaches zero at, the reflection level. Since more slowly 
moving waves spend a longer time in a given location, and are thus able 
to exchange energies more effectively, the decay of the pump wave 
into Langmuir ((n^) and ion acoustic (cd^) waves occurs mainly near the 
reflection level. Figure 3(u) shows two possible synchronism parallelo- 
grams denoted by subscripts 1 and 2, 

= “O, ■ ’ \l,2 = ^ 1,2 * *^0 ’ 

Where « k^, kp and may thus be approximated to vanish near the 
reflection level, and ^a^ satisfy the dispersion 

relations 



cn 


k v 
a 


a 


( 2 ) 


Here k^ is the electronic Debye wave-number: is the angular 

cyclotron frequency; © is the angle between and the geomagnetic 

field; and the acoustic speed is given by v = a[(3 0 + 

where a. is the square root of the ion/electron mass ratio, m is the 

electron mass and ^ are the electron and ion temperatures in energy 

units. 

At Arecibo, the intensity of this enhanced Langmuir wave is measured 
by sending a diagnostic wave vertically and determining its scattering 
intensity. Figure 3(^) illustrates the scattering synchronism conditions 
among the incident (tn^) , the scattered (down-shifted, cWjj , and 
up-shifted, tUy) , and the Langmuir waves, 

= “il,2 T n.1,2 ^ \u = HLi, 2 - \i,2 • (3) 


The dispersion relations for scattered and incident waves can be 
approximated by o) = kc , where c is the speed of light. 



CJ 



FIG, 3- Three-wave interactions of backscatter (dashed 
parallelograms) and heating (solid parallelograms). 

(a) The incident UHF waves and scatter off upward- 

propagating > O) and downward-propagating < O) 

Langmuir waves (l), at the same height, to give down-shifted 
(cWjj) up-shifted (co^) lines. 
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FIG. 3. Three-wave interactions of backscatter (dashed 
parallelograms) and heating (solid parallelograms). 

(b) The incident UHF wave (cn^) scatters off an upward- 
propagating Langmuir wave > O) to give a down- 

shifted line (cDj^), at a height above that at which it 
scatters off a downward -propagating Langmuir wave 
k^ < O) to give an up-shifted line . 

T 



If we 


Figures 1 and 2 are drawn against Af = foi* the 

down-shifted and Af = (f . + f.) - f for the up-shifted lines, 
neglect for a moment the effects of of f -synchronism and ionospheric 
inhomogeneity^ the scattered waves peak when and meet the 

synchronism conditions (l) and (3)^ i-e-j when Acq = m ^ . In order 

to obtain the value of q> ^ we first note that use of the dispersion 
relations of the scattered waves (co = kc) reduces ( 3 ) to 






^ Gl) 




W 


Since 2Jt X (5-6 MHz)] « [= ^ 2^ X (^30 MHz)] for the 

Arecibo experiment^ we may drop from (4). Therefore, . 

Assuming k^ 0 in (l), and using the ion acoustic wave dispersion 
relation in ( 2 ), reduces ( 4 ) to 


CO 


al 


a) 


b2 


2v 

a 

CO. 

c 1 


(5) 


For the temperature ^ 1200^ K and a I/I 7 O, the scattered 

waves then peak at Af = ^al ^ ^a2 ^ kHz which agrees with Figure 1. 

We now consider the effect of the second term on the RHS of (4). 
Due to the second term, f . The frequencies of the two 

Langmuir waves, on the other hand, are almost the same, ss co^ = co^ 

from ( 1 ) because . Therefore, the up-shifted wave must be 

scattered at a height different from that for the down -shifted wave. 
This is illustrated in Figure 3(b) . These two scattering heights can 
be obtained from the Langmuir wave dispersion relation ( 2 ) by assuming 
a linear electron density profile, ^ geomagnetic 

field, and uniform electron temperature, and by substituting from ( 4 ) 
to obtain 



12(n^(cn^ T cn^) 




( 6 ) 
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Here, z is the height measured from the base of the ionosphere, h 

1/2 ^ 

is the scale height, and v [ = (© /m) ] is the electron thermal velocity. 

© © 

In the last -term, cup = is taken because io « cn. . For typical 

■''to a 1 

values in the Arecibo experiment of f = 5*62 MHz, f =1.05 MHz, 

I O n 

0 = 40 , and h^ = 3 X 10 m , the up-shifted wave is scattered at 
about 900 m below the reflection height, and is at about 10 m above the 
height where the down-shifted wave is scattered. This height difference 
for the down and up-shifted waves may seem minor. As we shall see later, 
however, the intensity of the enhanced Langmuir wave is a sensitive 
function of height, and this small difference in height results in a 
considerable difference in the amplitudes of the down- and up-shifted 
spectra. 
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3. Scattered Spectrum Intensity 


The power scattered by the enhanced Langmuir waves is given by 
(Harker^ 1972) 


dP 

s 


0 T T ‘ 

(h + z)^ 



( 7 ) 


where and k ^ are the local electron density and the Debye wave- 

0 CIO 2 

number at the reflection level, T*9^ ^ n' ) is "the Thompson 

backscattering cross-section, P,p( = 2.5 is the diagnostic beam 

power, a(;= 150 ir m ) and ^0^) are the antenna aperture 

and efficiency, respectively, h(= I 50 km) is the height of the base of 

the ionosphere from the ground, and 1 is the steady state spectral 

jK 

energy density of the Langmuir wave, which can be obtained from the 
well known kinetic wave equation of parametric interaction among the 
pump (ooq), Langmuir (cn^), and ion acoustic (cd^) waves. Changing variable 
z to by (6), dz/h^ ~ 2a3^dLO^/LD^ , and substituting dP^/do 3 ^ - 

and (4) into ( 7 ); gives 


T 

s 


c^k ^ 
B dO 





h + z 



( 8 ) 


where is Boltzmann's constant. 

So far, we have considered the case where the synchronism condition 
( 1 ) is met precisely. We will now include the effect of a slight fre- 
quency mismatch , By neglecting higher order processes, 

such as nonlinear ion Landau damping or electron orbit perturbation, I 

^ k 

may then be written in a simple form (Barker, I 972 ) . 




PS0cot 0 


where P and the resonance function 0 are given by 


(9) 
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p = 


16 n^© 



0 = 


1 + (6m/Y )‘ 

a 


( 10 ) 


Here, y^{= 650/sec) and Yg^(= 8000 jr/sec) are the linear damping rates 
of the electron plasma and ion acoustic waves, respectively, and 
is the pump field intensity at z 0 , which is given by 


=0 = 




e^jtch 

0 


^ „2 

2 0 


130 , 


( 11 ) 


where is the HF pump power in watts, and G is the antenna gain. 

All we need to know now is the swelling factor, S , in ( 9 )* Then, 
the RHS of (8) can be evaluated. In what follows, we neglect the geo- 
magnetic field. After the pump field is found, however, it is regarded 
as being aligned in the direction of the geomagnetic field. In the 
absence of nonlinear absorption, the amplitude variation of a vertically 
propagating ordinary wave is customarily obtained by neglecting the 
Airy integral function, Bi , which grows with positive argument, because 
the field has to vanish as z "♦ «> . The field is expressed by the 
decaying Airy integral function, Ai ( Ginzburg, 1964 ; Budden, I96I) , 

In the present situation, where most of the pump wave energy is deposited 
at or near the reflection level, however, we assume that the electric 
field vanishes for z ^ h^ due to absorption, and phenomenologically 
keep the Bi term in the region below the reflection level (O ^ z ^ ^ 0 ^ ' 
In circuit terminology, this is analogous to a transmission line with a 
dissipative terminal impedance at z = h^ , whereas the usual treatment 
is analogous to an infinitely long nonuniform line. The solution to 
the wave equation in the ionosphere is then obtained by a similar 
procedure to that used in the absence of nonlinear absorption ( Ginzburg, 
1964 ; Budden, I96I) . We obtain 
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|e|^=s|Eq|^ , (o<z<hQ) , 



where the effective reflection coefficient, R , is measured at z - 0 
and C 2 nd Cq defined by 

The familiar results for S and R 
are recovered from (12) by letting the Bi term equal zero i^ uinzpurg ^ 

1964 ; Budden, I 96 I ) y as 

S = Al^(-C) , R = expj^-i|^ - ij] . (14) 

Our interests are in the scattering region. In this region 
1 6] » 1 and (3-2) becomes 

s = ^ |l + |r6 + 2|r|cos[1(c 3/2 _ j3/2j ^ ^ 

where is the phase of the reflection coefficient. The swelling 

factor is plotted in Figure 4 for |r[ =3 and \j; = 0 . 




=0 = (^f 


( 13 ) 


in the absence of the absorbing layer 


32 



M 

UO 



FIG, 4. Swelling factor for [ r] =3 and \[( = 0 as function of z/h^ , [The frequency 

scales are obtained from ( 6 )]. 


Now, the temperature spectrum (8) can be evaluated by substituting 
from (9) -(11) and (I 5 ) . Here, we note that I in (9) is a function 
of z primarily through y ^ y ^ ♦ Again changing variable 

z to by (6) , is a function of f only. The spectrum is 

plotted against ^f(^ f^ - f) in Figure 5 , with the phase of the 
reflection coefficient as a parameter. 
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TEMPERATURE (10^ '’K) 



FIG. 5* Backscatter spectra. (a) i|f =.-jt/2^ 0 . 
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FIG, Backscatter spectra, (b) i|r = jt/2 
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4. Asymmetry of the Lines 

The spectra in Figure 5(^) show an overall agreement with the 

experimental results of Figure 1, Each spectrum has peaks near 

Af - f = .4- kHz where the frequencies match precisely (f - f^ - f ) , 
a 0 a 

These peaks are due to electron plasma waves parametrically enhanced by 

the pump wave. The peak temperatures [ (6 13) X 10^ K] are in the 

range of experimental values [ Kantor^ 197^]* There are several maxima 

and minima near the peak of each spectrum. These are due to the swelling 

factor. The frequency intervals between adjacent maxima (10 l4- kHz) 

are in good agreement with the experimental results in the range 

8 18 kHz ( Kantorj 19T^j Figures 3 and 9)» The intensities of the 

up-shifted and down-shifted peaks differ from each other. This is due 

to the 10 m difference in heights (about 1/10 of the distance of the 

adjacent maximum -minimum of the swelling factor) where the scattered 

waves of f = f f are generated, as calculated by (6), For tjr = 0 , 

0 a 

Figure h shows positive slopes of S at both heights. Since the down- 
shifted plasma line is scattered from a height higher than that for the 
up-shifted, the down-shifted line strength is higher than the up-shifted * 
line, in agreement with the experimental results. 

For i|j = -it/2 , S shifts, and the heights for the down- and up- 
shifts are again on the positive slope somewhere between A and B in 
Figure 4 where the gradient of S is steeper. Therefore, the down- 
shifted line is stronger than the up-shifted again, but the difference in 
strength is greater than for ^[f = 0 . Figure 5(^) shows this difference 
as about 29 ^ of the up-shifted line amplitude, which is much smaller 
than that shown in Figure 1, However, it depends on the strength of 
the up-shifted line in a quite involved way and most of the experimental 
data indicate smaller differences than that shown in Figure 1, [See 
Kan tor (19/4, Figure 9) Carlson (1972, Figures 6 and 7)1* 

For \|r ^ jt/2 ' and it , on the other hand, the heights for the down- 
and up-shifted lines are on the negative slope somewhere between C 
and D in Figure 4. The down-shifted peaks are therefore lower in 
amplitude than the up-shifted, as shown in Figure 5(^) • 



5 . Frequency Displacements 


Since we expect the peak of the spectri4m to occur when I is a 
maximuTti_^ we can estimate the displacements of the plasma line frequencies 
from the pump frequency^ Af , from ( 9 ) . Now I is maximum when S0 
takes its maximum value. This occurs when d(S0)/da:»^ = 0 . Use of 
( 10 ) then gives 


— 

S doojj 


26cd 



( 16 ) 


We now expand S about the frequency precisely satisfying the synchro- 
nism condition (l), i.e., 6oci“0. 


dS^ , d^S 


dco^ 


To first order in ^ (16) then gives 


6(jd - 


0 

Sq 0^1 


1 


2 d 

Y ^ H ^ 


( 17 ) 


( 18 ) 


2 

Since is related to z through (6), d(z/h^)/dm^ - 200 ^/ 0 )^ « 2/(j0q ^ 

the frequency displacement of the Langmuir wave from the pump frequency^ 

Ao:) “ “ oop = tn - 6oi , becomes 

0 -L a ’ 


AjD = CD 

a 


1 ■ 2( Y^/o)q) ( Sq/Sq) 


( 19 ) 


where the prime represents a derivative with respect to z/h 

Now consider the dependence of Ao) on ij; from (I 9 ) • 

to know is the behavior of s'/S , Since S varies with 

^ ^ 6 —1/2 

through the cosine function in ( 15 )> fc>r j r| = 3 ^ 


0 * 

All we need 
z primarily 
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= -*'fW 


F(^) 


Sin ^ 

1/2 ^ 

2/3 + cos» 


(SO) 


- 3/2 3/2 

where = i|-(^ - Q )/3 + is the argument of the cosine function 

^ / 1/2 
in (15)^ and its ( z/h^) -derivative is ^ = ( 2 ojQh^/c)(l - z/h^) 

Since the variation of ^ ^ with respect to z is negligible, differen- 
tiating ( 20 ) once again gives 



/2 dF 


( 21 ) 


The function F is shown in Figure 6 in relation to , the origin 
being chosen so as to make ijf the abscissa. Except near the minima of 
Sq , which are not of interest because of the small pump amplitude, F 
is nearly linear with \j;, with slope about 1/2. For h^ - z 9 OO m, 

~ ijf 746 X 2it - 0.3 and therefore F may be approximated as 

F i (t - 0.3) , (for If I < 0.8jt) . ( 22 ) 

Since ^ 1224, ( 2 O) and (2l) reduce to 


- 6l2(>li - 0.3) , - 7-5 X ic/ . (23) 

Equation (I 9 ) therefore becomes 

Af f + 3.7 kHz . (24) 

We now see from (23) and (24) and from Figure 6 that, for -n < f s: 0.3, 

the slope of is positive and Af < 4 kHz. For 0.3 < f ^ Jt , on 

the other hand, has negative slope and Af > 4 kHz. 

Figure 7 shows Af from (24), and more accurate numerical values 

for f .= -jt/2 , 0 and it/2 obtained from Figure 5- These show that 
our results are in good agreement with Kantor’s observation 
Af = 3 4 kHz when -it/4 ^ ^ 0 , where the pump amplitude is nearly 

maximum and its slope is positive. It should also be noted that the 
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FIG, 6. Swelling factor_, , and Function F. 
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iiAf(kHz) 
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X DOWN-SHIFTED x 
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- tr -tt/2. 0 7t/2 


FIG, 7* Variation of frequency shift with phase 
of reflection coefficient. 
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frequency displacements of the down-shifted lines are always larger 
than those of the up-shifted, in agreement with the results of Kantor 
shown in Figure 2. The frequency difference is independent of ijr , 
with a value about O .5 kHz. Again, this agrees with the experimental 
value of 0.4 ± 0.1 kHz. 

The schematic diagram in Figure 8 offers a simple pictorial 
explanation for the larger displacements of the down-shifted lines than 
those of the up-shifted. We first note from (23) that S is convex 
throughout the region, except near its minima. Since slopes of S at 
the scattering heights are positive for -n < \|; ^0.3 and negative for 
0.3 < Jt , they may be drawn against Af as in Figure 8(b) and (c). 
Note that, since Af = f^ — f , the slopes of S in (a) are reversed 
in (b) and (c). The figure also shows the resonance function 0 
peaking at Af = f^ . Then, Af of the peaks of the product 0S 
increases or decreases from f^ , depending on the sign of the slope 
of S . The figure clearly demonstrates that, for both cases of (b) 
and (c), Af of the peak of 0S is larger for the down-shifted than 
for the up-shifted line. 
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FIG. 8. Frequency displacements of plasma lines 

[Schematic, Four values of swelling factor indicated 
in (a) are all at Af = f^ , In (b) and (c)_, swelling 
factors for up-shifted (s^) and down-shifted (S^^) 
lines near Af = f^ are approximated by straight lines.]. 


23 


6 . Conclusion 

Since the scattered spectrum intensity is proportional to the pump 
field amplitude, the location where scattering occurs in relation to 
the standing wave pattern of the pump significantly affects the spectrum. 
Because of the random nature of the ionosphere, the peak amplitudes of 
the plasma lines and their frequency shifts vary widely with time. 
However, their average features agree well with the present theory if 
scattering occurs where the pump field is nearly maximum and its slope 
is positive. This happens when the phase of the reflection coefficient 
is confined within the range ijf ^ -Jt /4 0 . It is therefore a subject 

of further experimental verification to establish whether ijf averaged 
over the time period of an observation is indeed near or within this 
range. 

The several peaks in the region near the HF transmitter frequency 
have been explained previously as the transfer of power from one 
Langmuir wave to another at a slightly lower frequency via an ion 
acoustic beat wave ( Kuq and Fejer, 1972, Kruer and Valeo, 1973; Fejer 
and Kuo, I973) . The present analysis shows, however, that the standing 
wave pattern of the pump wave can also result in these pumps. Further 
work is required to determine whether the peaks actually observed are 
due to the combined effects of the two mechanisms, or predominantly due 
to one of them. 

Acknowledgement . This work was supported by the National Aeronautics 
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